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Abstract 



We show that the distribution of the coefficients of the g-derangement numbers is 
asymptotically normal. We also show that this property holds for the g-derangement 
. numbers of type B. 
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(Si ' 1 Introduction 

ON ' 
O 

Let ©n denote the symmetric group of permutations on [n] = {1,2, ... ,n}. Let @> n denote 
\Q ■ the set of derangements, i.e., 

o 

OO 

O 



@n = {tT = 7ri7T2 • • • TT n £ & n [ 7Tj / i, % = 1, 2, . . . , Tl}. 

The major index of a permutation ir = iri^ ■ ■ ■ vr n is defined by 



^ 1 maj(7r) = ^ i. 

7Tj>7Ti+l 

The following formula was derived by Gessel and published in |llj : 

dn{q) = £ ^ ajW = W E TfT^' M 

7re^„ fc=0 ^ 

where [0], = [0],! = 1 and for k > 1, [fc], = l + g + g 2 + - • and [k) q \ = [k] q [k-l] q ■ ■ ■ [l) q . 

The coefficients of d n (q) are given in Table [TTTI for n < 6. Combinatorial proofs of (jl.ip have 
been found by Wachs [IB], and Chen and Xu [BJ. 

In this paper, we will show that the limiting distribution of the coefficients of d n (q), that 
is, the major index of a random derangement, is normal, see Figure 11.11 Moreover, we will 
show that the limiting distribution of the g-derangement numbers of type B is also normal, 
see Figure 11.21 
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Table 1.1: The g-derangement numbers of type A for n < 6. 




Figure 1.1: The distribution of the coefficients of @xo compared with the normal distribution. 



Write the set {1,2, ... ,n} as [n\. Let 6^ denote the hyperoctahedral group of permu- 
tations on [n] U [n], called signed permutations or 5 n -permutations, see Bjorner and Brenti 
[5]. Let f^f denote the set of B n - derangements on [n], namely, 

@n = = ^1*2 • • • 7T„ € : TTj ^ t, » = 1, 2, . . . , n}. 

For example, = {1}, &2 = {^2, 21, 21, 21, 21}. For i? n -permutations, Adin and Roichman 
[3] introduced the notion of the flag major index, or the fmaj index for short, defined by 

fmaj(7r) = 2maj(-7r) + neg(7r), 

where maj(7r) is the major index of ir with respect to the following order on [n] U [n]: 

n<---<2<I<l<2<---<n, 

and neg(-7r) is the number of 7Tj's in [n], see also Adin, Brenti and Roichman [2], and Chow and 
Gessel [8]. For example, the flag major of the Z^-permutation 3512674 equals 2 x 11 + 3 = 25. 
Chow [7] derived the following formula for the g-derangement numbers of type B: 

d B M = E ^ fmajW = [^» E ( U „ , (i.2) 
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Table 1.2: The g-derangement numbers of type B for n < 4. 



where [2&] g !! = [2/c]q[2/c — 2] g ■ ■ • [2] g . For n < 4, the coefficients of the polynomials d^(q) are 
given in Table [L2l 

Based on the formula (|1.2j) , we will show that the limiting distribution of the coefficients 
of d^(q) is normal. Figure [l~2l is an illustration of the distribution for n = 10. 




Figure 1.2: The distribution of the coefficients of compared with the normal distribution. 



2 The Limiting Distribution of the Coefficients of d n (q) 



The aim of this section is to show that the limiting distribution of the coefficients of d n (q) is 
normal. We write 

1, if k = n; 



[n] q [n - l} g ■ ■ ■ [k + l] q , else. 



/n,fe(?) = < 

Then we can express d n (q) as 

n 

dn(q)=Y,(- 1 ) k( l Q fnAQ)- 



(2.1) 



k=0 



Let D n = \3> n \ be the number of derangements in Q} n . For example, D\ = 0, D2 = 1, -D3 = 2, 
£> 4 = 9, D 5 = 44. 



We will adopt the common notation in asymptotic analysis. If f(n) and g(n) are two 
functions of n, then 

• f(n) = o(g(n)) means that lim n ^oo f(n)/g(n) = 0; 

• f( n ) ~ g(n) means that lim^oo \ f(n)\/\g(n)\ = 1. 



We now recall some basic facts about the derangement numbers D n , see, for example, 
Stanley p2]. For n > 3, 



n 



E 



(-i)* 



k=0 

(n -l)(£>„_i + D n _ 2 ) 
l 



n! 1 
7 + 2 



n! 



(2.2) 

(2.3) 
(2.4) 

(2.5) 



where the symbol [a^J denotes the largest integer not exceeding x. From (|2.3p it immediately 
follows that 

D n -i 1 (-l) n 1 , v 

+ o(l). (2.6) 



D n n nD n n 



While it is common to use maj(7r) to denote the major index of a permutation tt, there 
does not seem to be any confusion if we also use maj to denote the major index of a random 
derangement on [n]. The probability generating function of maj is clearly d n (x)/D n , whereas 
the moment generating function of maj is given by 



M n {x) = ^d n {e x ) = ±-[n] e AY j 



^ (-l) fc e*(5) 



fc=0 



[fc]e«! 



(2.7) 



1/2 

Let E n , V n and a n = V n denote the expectation, the variance and the standard deviation of 
maj respectively. Then the probability generating function d n {q) of the normalized random 
variable (maj — E n )/a n equals 

dn{q) = Y q^'^~ E ^/ Un =q- En / an d n {q l ^ n ). 

Thus by the definition ()2.7|) . the moment generating function of (maj — E n )/a n equals 

Mn{t) = dnie^/Dn = exp(-t E n /a n )M n (t/a n ). (2.8) 



2.1 The expectation and variance 



We now compute the expectation and variance of the major index maj of a random derange- 
ment on [n]. 
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Theorem 2.1 The expectation E n and variance V n of the random variable maj given by 

lfn\[ D n _ 2 \ n 2 -n + l (-l) n (n - 1) , x 

and 



Vn 



2ra 3 + 3ra 2 - 5n - 16 9n 3 - An 2 - 46n + 41 (-1)" (n-l 



72 



144 



A, 



4D« 



(2.10) 



Here we give only a sketch of the proof, and detailed steps are omitted. 
Proof. The generating function (jl.ip implies that 



Dr. 



Da ' 
2 <(1) 



(2.11) 



(2.12) 



where <i^(g) and d'^(q) are the first and second derivatives of d n {q). From (12. 1|) and (12, 2p . we 
find 

1 



^(l)=E(-l) fe [(2)/(l) + /'(!) 



2 V2 



(£> n + L» n _ 2 ). 



So ([13]) follows from ([2TTT]) . Differentiating (|2TT|) twice yields 



= [(*) ((*) - 1) m +2Q) m + rci) 

The following relations can be easily verified: 

n 



(2.13) 



fc=0 



3 /n 



2 V3 



(n + l)D n _ 2j 



k=0 ^ 3 ' 



Now, using (|2.13p and (|2.4p . we deduce that 



1 In 



According to (|2.12p , 
n 



144 



-(n-l) (27n 2 - 13n - 23) -£=i + (9ra 3 + 13n 2 - 7n - 38) 



n 2 n(n — 1) -D n -i 

T I aT 

In view of fl23D and (pUj) . we obtain OTTO]) . 
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We note that the formula (|2,9p for the expectation of the major index can also be derived 
by a combinatorial argument, the details are omitted. Based on the estimates (12, 5j) and (12. 6ft . 
we derive the following approximations. 



Corollary 2.2 We have the following asymptotic estimates: 



n 



- n 1 , . n 3 n 2 5n 2 



^ = T"4 + 4 +0(1) ' K= 36 + 24 - 72 - 9 +0(1) - 
2.2 The limiting distribution 

It is well-known that the moment generating function of a random variable determines its 
distribution by Curtiss's theorem (see Curtiss [9] or Sachkov [H]). In particular, if the 
moment generating function M n {x) of a random variable £ n has the limit 

lim M n (x) = e x2/2 , 

then £ n has as an asymptotically standard normal distribution as n trending to infinity. 

We will need Tannery's theorem (see Tannery [13]) which is essential in the proofs of 
Lemma 12.41 and Lemma [32 



Theorem 2.3 (Tannery's theorem) Let {ffc(?i)}fc>o be an infinite series satisfying the fol- 
lowing two conditions. 

• For any fixed k, there holds linin^oo Ujt(n) = w k - 

• For any non-negative integer k, \vk{n)\ < M k , where M k independent of n and the 
series Ylk>o is convergent. 

Then 

m(n) oo 

lim Y] v k (n) = y^ w k , 

n— >oo * — * * — • 

k=0 k=0 

where m(n) is an increasing integer-valued function which trends steadily to infinity as n 
does. 



Lemma 2.4 For any \x\ < 1 and bounded \t\ < M, we have 

n u 
X 



lim ErZi \ = e% - ( 2 - 14 ) 

Proof. We apply Tannery's theorem and set 



k 

X 



v k{n) = rn 

[k\ e -t/a n \ 
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and m(n) = n. Then for any fixed fc, by Corollary 12. 2[ it is clear that 

x^ 

w k = lim v k (n) = —. 
Note that the right hand side of (|2.14j) can be expressed as 

OO oo 



OO 1, oo 
x ^ " 

' k\ 

k=0 k=0 



By virtue of Tannery's theorem, to prove (|2.14p it suffices to find an upper bound M}~ for 

| X^ | 

\vk(n)\ 



[k] e -t/cr n 



such that Mfc is independent of n and X^£=o converges. We claim that there exists a 
constant c G (0,1] such that M*. = (1 + c) 1_fc is the desired upper bound and this bound 
clearly implies the convergence of Y^=o = l/c + 2 + c. 

For t < 0, we have e _ '/ cr " > 1 and thus 



\x k \ \x k \ I I 



< 



< — < —r— < M k . 

U — h — 9fc-l - fe 



[fc] e -*/«rn! " k\ ~ k\ ~ 2 k 

For t > 0, Corollary 12.21 implies that <r n has a positive lower bound as n runs over all positive 
integers and so does e _t// °" n . Suppose that e~ t l (Jn > ct £ (0, 1]. Since the function e _t// °" n is 
continuous in t and t is bounded, there exists a constant c £ (0, 1] independent of t so that 
e -t/f„ > c £ or |^| < ^ Hence for any k > 1, 

■ jU I fc II 

kc tt kc 



[ fe ] e -*/*n ! ^ 1 + e -t/<7„ _| |_ e -U-l)t/<r» 



< TT — = 



This completes the proof. I 

In the computation of the moment generating function of maj , we will need the Bernoulli 
numbers B k which have the following generating function, 



k=0 

The first few Bernoulli numbers are 



X \ X , 

— = E^i[- ( 2 - 15 ) 



B = 1, B x = -1/2, B 2 = 1/6, 5 3 = 0, B A = -1/30. 
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Moreover, -621+1 = for any i > 1. Alzer [3] establishes sharp bounds for \B% n \ leading to 
the following asymptotic formula (see also [U pp. 805]) which will be needed in the proof of 
Lemma 12.51 

2- (2n)! 
(2ir) 2n 



1^1 ~ ( 2 - 16 ) 



Lemma 2.5 For any bounded \t\ < M, we /love 



i im y f^t y , 2i _ v = 



where B<n are the Bernoulli numbers. 



Proof. Let a, /3 and 7 be three constants such that a > 1, /? > 36, and < 7 < 1/2. Let ./V 
be a fixed integer satisfying the following three conditions: 

• n + 1 < an for any n > N; 

• er 2 — n 3 /(3 > for any n> N; 

• 27ri\H/ 2 > Mav73. 



The existence of such N is obvious. Let i > 2 and n > N. From the inequalities 



n r n+l 

E 2i - 1) < / (* 2i - 1 
.2 ^ ^3, 



n+1 , w N , (n + l) 2i+1 - 1 (n + l) 2i+1 (cm) 2i+1 

1 Mi = n < < 

2i + 1 5 5 



and the assumption <r 2 > n //?, we deduce that 

/3< M"* 1 _ a (a^'f 



|r£(i 2i -i)< 

In light of the inequality 



(7~* ' — ' N n 3 * 5 5 n 1 1 



1 1 

< 



n7 i n (l- 7 )i-l j\r 7 i n i-2 7 ' 
we see that 

i im v ^ 2i * 2t v d 21 -i\<-]im ( v 1^1 ("v^y A 2 7 -i f2 18) 

^ (2i ) (2i)! a 2 / ^ [J ' ~ 5 ™ ^ (2t) (2i)! N* J ' [ ' 

By the asymptotic estimate (|2.16p for Bernoulli numbers, we see that the radius of conver- 
gence (see, for example, Howie [10]) of the series on the right hand of (|2.18|) equals 

lim ( \m b^_ V* = ^ >M . 

i-oo^(2t)(2»)! J ay/ft 

Since lim n _ +00 n 27_1 = 0, we conclude that the series in (|2.17p is absolutely convergent to 
zero for \t\ < M. I 
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The following lemma gives an expression of the moment generating function of the 
random variable maj in term of the Bernoulli numbers. This lemma will be needed in the 
proof of Theorem 13.61 



Lemma 2.6 The moment generating function of maj equals 



M n {x) - n exp | 4 +^ (2 . )(2j)! ^ U lj I ^ 



Proof. By the formula (12. 7h . we need to express [n] e »! and e x ( 2 ) /[/c] e *! in terms of Bernoulli 
numbers. It is known that, see, for example, Mcintosh |12j . 

l_ e =x-exp 



fc-Jfc! 
fc=i 



Thus for any j > 1, 



1-e^-xj-expH — ^ =-xj-exp — + ^ 



b]ex = =j-^[ —r~ + g ( 2.) (2.)! J • 



Therefore, 



is*! 



lJ>]e» = 11 J- exp + ^ 



(2i) (2i) 

3=1 j=l \ i=l y 1 K ' 



n\- exp 



i=l v ' v 7 j=i 



Observe that 



„m / A 1 - \ A 1 - e x e JX 1 



n i-e^ ) n 1 _ e ^ e x \kh- x \' ^ 2 ' 20 ^ 



[fc]e*! I ■ IJj l-e x J 1 _ e jx e x [ A ] £ 

Substituting ()2.19p and (12.200 into (j2.7l) . we obtain the desired expression. 



Theorem 2.7 Zei maj fee the major index of a random derangement on [n]. Then the dis- 
tribution of the random variable 

maj - E n 

u = 

converges to the standard normal distribution as n — > oo. 
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Proof. By Curtiss's theorem and (|2.8p . the normality of the distribution of the standardized 
random variable £ n can be justified by the following relation 

lim e- tEn / an M n {t/a n ) = e* 2 / 2 . 

n^oo 

By virtue of Lemma 12.61 the above relation can be restated as 

lim nL exp ( J3l + + y B * t7A y d 2t { ~ 1)k - e t2 ' 2 

■SSo D n 6XP y a n + 4a„ + ^ (2i) (2i)!a* ^ V ^ J ^ [fc]^ ! " ' 

First of all, the estimate (|2.5|) implies that 

lim n!/D n = e. (2.21) 

n— »oo 

By Corollary 12.21 for bounded i we have 

fn(n-l)t tE n \ , i /n(n-l) \ 
lim -i- i = lim — -i— - -E n ) = 0, (2.22) 



It is easily checked that 

1 n 

lim — V (j 2 - 1) = 12. 
n— >oo rrs *- ■" ' 



<7„ . 



In view of Lemma 12.51 and the fact that B2 = 1/6, we have 

lim y f** _ y (f -!)= l im y (f - l) = L. (2.23) 

i=l K ' x ' n j=l " j=l 



Finally, taking x = — 1 in Lemma 12.41 we get 

■^Spfcb- 8 "'- (2 ' 24) 

Combining (I23TD . (12321) . (12331) and <[233j) . we complete the proof. I 

3 The Limiting Distribution of the Coefficients of d^(q) 

In this section, we show that the limiting distribution of the ^-derangement numbers is 

normal. Let be the number of -B n -derangements on [n]. The first few values of are 

Df = 1, D% = 5, D§ = 29, Df = 233, D§ = 2329, D§ = 27949. 

For n > 3, we have 

n C— Fl fe 

Ztf = c£(l) = (2n)!l£^L (3 .1) 

= 2nZ#_ 1 + (-1)" (3.2) 

= (2n - + (2n - 2)£>£_ 2 (3.3) 



(2n)l\ 1 
- — 1 — 

v 7 ^ 2. 



(2w)!! 



(3.4) 
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For completeness, we present a proof for ([3 



u— n \ / ' \ i — „ i 



> , (2fc)!! 

fc=0 v ' \ k=n+l v ; 



It is easy to see that the absolute value of the remainder 

-If 



(2Jfc)!! 

k=n+l y ' 

is not greater than the absolute value of the (n + l)-st term of the alternating series, i.e., 
l/(2n + 2)!!. This yields 

D* = (2n)ll/Vi-y n , 

where 

\y n \ = \{2n)\\r n \ < (2n + < 1/4. 
Since is an integer, (|3.4|) is verified. 



+ o(l). (3.5) 



From (|3T2j) it follows that 



D§ 2n 2nD§ 2 



n 



Let , Vn and = (V^) 1 ^ 2 denote the expectation, the variance and the standard 
deviation of fmaj respectively. We also use fmaj to denote the fmaj index of a random 
-B n -derangements on [n]. The probability generating function of fmaj is 

The moment generating function of fmaj is given by 

(*) = 7^^n (el = E ( 7^ E fma j^)" J • ( 3 - 6 ) 

The normalized random variable (fmaj — E^)/a^ equals 

MB(t) = exp [-tE*/a*) M* (t/a*) . (3.7) 

3.1 The expectation and variance 

Let [x l ]f(x) to denote the coefficient of x % in the expansion of f(x). Then the expectation 
and variance of fmaj can be expressed in terms of the moment generating function M,^(x): 

^ = ^ E fma J = \*\M%[x). (3.8) 

V n = ( 4 E fm& j 2 J " = 2 M M n (*) ~ i^f ■ (3.9) 
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Let {x 2 ) f(x) to denote the truncated sum of f(x) by keeping the terms up to x 2 . Once 
{x 2 ) M^(x) is computed, then the first and the second moments are easily extracted. In this 
notation, we have 

( x 2 ) e rx = 1 + rx + r 2 x 2 /2. 

Moreover, 

2j-i 



r=0 ^ ' 



j=fc+l r=0 i=fc+l 



where 



ci 



(2n)!! / n 2 -*; 2 2 



j(2j - l)(4j - 1) 1 , (2i\ (2j\ £ J_ 



2^ 



6 2 j ^ \2 \2 2i2j 

j=k+l k+l<i<j<n v 7 x ' 

_ (n - A;)(9n 3 + 4n 2 + 9fcn 2 + 6n - 9fc 2 n + 4/cn - 1 + 6/c - 9& 3 + 4fc 

~ 72 ' 
By the definition (|3.6p . we find 

n r , .. n n 2j-l 

" fc=0 1 Je:E " n k=0 j=k+l r=0 

It follows that 

<* 2 >^ iMb A , ^ ^ , ^(fc-l) 2 2 \ (2n)!! / , n 2 -fc 2 



1 A, si .(2n)!! / (n 2 -k 2 ,, A „\ 



fc=0 

where 



n 2 -k 2 k 2 (k-l) 2 
C2 = Cl + k (k-l) + 



2 2 

Let (k)i = k(k — 1) • • • (A; — i + 1) be the lower factorial. We get 

9(k) 4 + U(k) 3 + (18n 2 - 27){k) 2 - I8n 2 k + (9n 4 + 4n 3 + 6n 2 - n) 



(>2 



72 

Combining (pU]l . and (pl8]l . we find 



n 2 n 



1 

n 


n 

£<- 

fc=0 


nfe (2n)H 
J (2A;)!! 


n 2 


3ro 

+ T 






) DB ■ 



+ k(k-l) 
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Now, the variance of fmaj equals 



fc T 9 14 18n 2 - 27 



D-d 



fc=0 



+ T^T^ "rrrr + 



2 4 (2/fc-8)!! 2 3 (2fc-6)!! 2 2 (2A:-4)!! 



18n 2 9n 4 + 4n 3 + 6n 2 — n 
+ 



2(2fc-2)!! (2fc)!! 



It can be deduced that 



2i iirB/^ _ n{72n z + 140ra 2 - 22n - 101) ra(216n 3 - 356n 2 - 186n + 127) D^_ ± 

576 576 W~' 



Theorem 3.1 The expectation E% and variance o/fmaj given by 



n 2 4 \ 2 4 J D% 
and 

B _ n (68n 2 - 40n - 101) n (72n 3 - 212n 2 - 78n + 127) n 2 (2n - 3) 2 / ' 



2 



288 288 D% 16 I 



In view of (|3.4|) and (|3.5|) . we obtain the following estimates. 

Corollary 3.2 We have the following asymptotic estimates: 

_, B n 2 3 d re 3 n 2 n 13 

E„ = — i h o(i), Vr = — i h o(i). 

2 8 v; 9 6 36 36 w 

3.2 The limiting distribution 

We aim to show that the limiting distribution of fmaj is normal. The following formula is 
analogous to Lemma 12.41 



Lemma 3.3 For any real x satisfying \x\ < 1 and bounded \t\ < M , 

n 

lim V - = e x ' 2 

n ^, n [2k] tlaB \\e kt l°* 

Proof. By virtue of Tannery's theorem, if suffices to find an upper bound M/% for 

I 7 1 ^ I 

\ v k(n)\ 



[2*] t/a BUe**/«tf 
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such that Mk is independent of n and Ylh=o converges. 

If t < 0, Corollary 13.21 implies that has a positive lower bound as n runs over all 
positive integers and so does e*/ ™ . Suppose that e*/ ™ > c\ £ (0, 1] for all |i| < M, where c\ 
is independent of t. Then for any k > 0, 



[2fc] e _ t/CT s!! e**/^ 2 (! + e -t/«£ + e -2t/o£ + . . . + e -(2j-i)t/<tf) e t/«r£ 



•J- 



= II fr < T » 4- I + ( -' fi -'-'./-^'. ) 

< (l + ci)- fc . 

Clearly, X^x)^ + °i)~ k * s convergent. 

We now assume that t > 0. Suppose e -i /°" n > C2 G (0, 1] where C2 is independent of t. 
Then for any k > 1, 

i ft ii 



_ = TT !±J 

[2fc] e _ t/CT B !! e fci /< 11 ^ + e -t/o£ + 1_ e -(5y-i)t/«£) e fct/wj? 

<II tt r 1 - 2Fi <a+^r fc - 

7_ n 1 + Co H h (V 



j=1 ± + c 2 H h c 2 

Similarly, ^Zjfc>o(l + c 2)~ fc is convergent. I 

The following formula, which is similar to Lemma 12.51 w iH be crucial in the proof of 
main theorem of this section. 

Lemma 3.4 For any bounded \t\ < M, 

lim y ^1 y f(2j) x - 1) = 0, (3.10) 

where B<n are the Bernoulli numbers. 

Proof. Let a, /5 and 7 be three constants such that a > 2, (3 > 9, and < 7 < 1/2. Let iV 
be a fixed integer satisfying the following three conditions: 



In + 2 < an for any n > N; 
(o-n) 2 - n 3 /{3 > for any n> N; 
2vriV7/2 > Ma^/fi. 
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The existence of such N is evident. Let i > 2 and n > N. We will show that the series in 
(|3.10p is convergent to zero absolutely. It is easy to derive the following upper bound: 

n n rn+ i 

((2j) 2i - 1) = 2 2i f l ~n<2 21 / t 2i dt 

3=1 3=1 Jl 

2 2l (n + l) 2i+1 (2n + 2) 2i+1 (an) 2i+1 

< — < < - — . 

2i + 1 5 5 



The rest of the proof is similar to that of Lemma 12.5 



Lemma 3.5 The following relation holds: 



M » w _ <^ - + g g ((2 ,, _ , g _JzlL. 



Proof. From (J3T6J) and (jl.2p . we have 



M n B (x) = d* {?) /D» = -g [2n] e *U { ~ * „ ■ (3.H) 



k=0 

Moreover, 



[2n] e *H 



TTF9-1 TTf9 -1 ^(2,-1) ^ £^((2^-1) 



i=i j=i 



and 

e xk 2 ^ * e (2j-l)x ^ * ^jx l_ e x l _ e -x j 

[2jfe] e *!! ~ JLA [2j] e * ~ JLJL ~ i _ e 2i^ ~ JUL 1 - e -2jx ~ [2k] e - x \\ ' (3 ' 13 ^ 

Substituting (j3. 12[) and (13.130 into (|3.1ip . we deduce the desired relation. I 



Theorem 3.6 Let fmaj be the flag major index of a random B n - derangement. Then the 
distribution of the random variable 

B fmaj - E% 
«n B 

n 

converges to the standard normal distribution as n — > oo. 
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Proof. By Curtiss's theorem and (|3,7p . normality of the distribution of fmaj follows from the 
identity 

lim e~ tE n/< M * (t/o*) = e t2 ' 2 . (3.14) 

By Lemma 13.51 the left hand side of (|3.14p can be expressed as the limit of the following 
expression: 

First, the estimate (j3.4|) implies that 

lim {2n)\\/D^ = y/e. (3.15) 

n— >oo 

By Corollary 13.21 for bounded £ we have 

/ tE% tn 2 \ , i /ra 2 R \ 
lim -- + = lim hr-^n =0. (3.16) 



»— x \ q-B 2a^ J n -> oo cr^ \ 2 



It is easy to check that 



((2i) 2 -l) = 12. 



hm 

n^oo ( rtB\ A 

\ u n ) j=l 

Based on Lemma 13.41 and the fact that B2 = 1/6, we see that 



i h, I ' '--■'>■ • B? 1 

3=1 

Finally, taking x = — 1 in Lemma 13.31 we get 



lim y ^il y ( (2 j) 2 * - 1) = hm ^ ? V ((2j) 2 - 1) = -. (3.17) 

"-o 2 . 2! (a*) 2 ^ 2 



lim y r n ( 1)fc - B = e- 1 ' 2 . (3.18) 

k=0 L J c 



-t/a. 



Combining (I3TT51) . (jgHgj) , (l3Tfl) and (pTTSjl . we obtain pTIj) . This completes the proof. I 
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